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Abstract 



We consider a random surface representation of the three-dimensional Ising model. The 
model exhibit scaling behaviour and a new critical index k which relates jstring for 
the bosonic string to the exponent a of the specific heat of the 3d Ising model is 
introduced. We try to determine k by numerical simulations. 
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1 Introduction 



The 2d Ising model formulated on the regular square lattice has a well know random 
walk representation. In fact we can write: 

Z(J3) = £exp(/3 J2 a n a m ) ~ £*(P) exp(-/i(/3)L(P)), (1.1) 

At(^) = -lntanh/3, $(P) = (-1)™ (P) . (1.2) 

In these formulas P is a closed (not necessarily connected) path without back- 
trackingQ on the lattice, L(P) the length of the path (i.e. the number of links 
constituting the path) and n(P) the number of times the path P self-intersects. 
The sum of paths exponentiate as a sum over connected closed paths and the free 
energy of the 2d Ising model: 

P(/?)~lnZ(/3)~ £ $(P)J_exp(-/i(/3)L(P)). (1.3) 

connected P \ J 

It is known that the 2d Ising model has a fermionic representation, a fact closely 
related to the sign factor $(P), and that the continuum field theory associated with 
the second order transition of the 2d Ising model is that of a free Majorana spinor. 

If we consider the 3d Ising model defined on a regular cubic lattice it has a random 
surface representation [H which can been viewed as the analogue of the random walk 



representation (|Ll)- (|1.2f) . It can be written as 



Z{(3) = £ exp(/3 J2 Wm) ~ E HS) ^ exp (-^)A(S)) . (1.4) 

CT (nm) S \ / 

//(/?) = -lntanh/3, ${S) = (-1)' (5) , (1.5) 

where denotes a closed lattice surface (without back-tracking) built from plaque- 
ttes, C(S) a symmetry factor for the surface @, A(S) the lattice area of S (i.e. 
the number of plaquettes constituting the surface S) while l(S) is the number of 
links where the surface intersects itself. (— l)"- 5 - 1 is obviously a surface analogy of 
the Kac-Ward factor (— l) n ( p ) for the random walk. Again the sum over surfaces 
exponentiates into a sum over connected closed surfaces @: 

P(/3)~lnZ(/?)~ J2 HS)-^-exp(-^)A(S)). (1.6) 

connected S v ' 



1 Back-tracking means in this context that two successive steps in the lattice path shares the 
same link. For surfaces the generalization is that two neighbouring plaquettes do no lie on top of 
each other. 
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Although it is by no means clear how to take the continuum limit of this lattice 
surface theory, due to the oscillating sign factor $(5), one would nevertheless think 
that the second order transition of the 3d Ising model should result in continuum 
field theory, and that this continuum theory should allow a representation like ( |1.6|) , 
only with continuum surfaces instead of lattice surfaces. For a closed path P in two 
dimensions it is possible to obtain a simple geometrical expression of the sign factor 
$(P): 

^(P)=exp(-i<fe i de t ), (1.7) 



.2 

where e*(£) is the tangent vector on the path P, but an analogous expression is not 
known for <&(S) for a surface S immersed in R 3 . One can nevertheless start on the 
lattice, introduce a fermionic representation of $(£) and take the formal continuum 
limit of the lattice fermionic action || || (see also || for related considerations). In 
this way the continuum limit becomes a string theory where the action is sum of an 
area term and the induced Dirac action plus an anomaly term on the world surface 
immersed in R 3 and we can formally view <&(S) as a result of first integrating over the 
fermions. If we accept this point of view, it makes sense to talk about a topological 
expansion of the surface theory, which at first glance might appear somewhat sur- 
prising since the original high temperature expansion over lattice surfaces involves 
crucial cancellations between surfaces of different topologies and cancellations be- 
tween orientable and non-orientable surfaces. If we assume that such a topological 
expansion is possible close to the critical point of the Ising modelQ we can discuss 
scaling properties of $>(S) which are consistent with this assumption and the known 
critical behaviour of the Ising model. 

Since the 3d Ising model has a second order phase transition at a certain value 
(3 C we have 

F(/3) = r.t(/3) + c(& - (3f~ a + Ls(/3) (1.8) 

where r.t(/3) stands for regular terms, l.s(/3) stands for less singular terms, while the 
critical exponent a is known from high temperature expansions and renormalization 



group study to be a ~ 0.11 or 0.125 |20, ETJ. How can such an expression match the 



left hand side of (|1.6| )? If we ignored sign factor <&(S) we would have an ordinary 

2 It should be stressed that it seems difficult not to have a topological expansion, possibly 
involving both orientable and non-orientable surfaces, if there really exists a continuum fermionic 
surface representation at the critical point. Viewed in this way the assumption of a topological 
expansion is more or less equivalent with the assumption of a continuum surface representation. 
We feel there are good reasons to expect the existence of such a continuum surface representation 
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bosonic string theory^. The entropy factor of bosonic surfaces of a fixed topology 
has (at least for c < 1) the form 

Z(A) ~ A^ 3 e^° A . (1.9) 
This leads to the well known singularity structure of the bosonic string: 

Z(fj) = £ Z(A)e-» A ~ r.t(//) + c(/i - /i ) 2 " 7 + Ls(/i) (1.10) 

which seems unrelated to ( |1.8|) , since /i > /x(/3 c ) . It is however intuitively clear 
that ($(5)) is small and it is reasonable to assume it has the following asymptotic 
behaviour as a function of the area A: 

, iE M^^.^ (1 .n) 

E{S\A(S)=A}e >loA{ - S > 

Under this assumption the rhs of (|1.6|) can be written as 

Z(A) $a e-» A ~ r.t( M ) + c(/i - (/x - /ii)) 2 " 7 " K + l.s(/x) (1.12) 

and comparing with fll.gj ) we reach the conclusion that 

« = 7 + k, (x - /jii = fx(p c ). (1.13) 



At this point it should be made clear that the same arguments can be used in 
the case of the ordinary bosonic random walk, $(P) and the 2d Ising model. It can 
be shown that we indeed have a relation similar to (|1 . 1 1|) : 

The exact values of fi and /2i depends on the discretized models used for the bosonic 
random walk. The result ( |1.14| ) is implicitly in the formalism developed for dis- 
cretized fermionic walks in ||, where the phase factor ( |1.7|) is represented as a 
product of Hilbert-Schmidt operators on S 1 . The largest eigenvalue of that oper- 
ator dominates in the scaling limit and one can prove that k = — 1 in agreement 
with ( p..l3| ) since 7 = 1 for a closed random walk in two dimensions and a = 
(logarithmic divergence) for the 2d Ising model. 

3 In the following wc assume the world sheet variables are somehow regularized, either by the 
original cubic lattice or by dynamical triangulations. 
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In the case of random surfaces a formula like ( |1.13| ) immediately leads to the 
interesting question of topology dependence: /i is independent of topology and it 
is reasonable to assume that the same is true for [ix- This assumption is needed for 
Ho — Hi = h(Pc) to make sense, since /3 C has no a priori relation with the topology 
We know on the other hand that 7 depends on topology in a linear way. For central 
charge less than or equal to one we have as a function of genus g of the surface: 

lg = lQ + K{c)g, K(c) > 0. (1.15) 

K(c) is not known for c > 1. In case the surfaces for c > 1 are dominated by so- 
called branched polymers we know that K(c) < for sufficiently large c |p2|, |23|. If 
K(c) < it is plausible that surfaces of spherical topology will dominate the critical 
behaviour, but this dominance could in principle be eliminated by a g-dependence of 
k. If that is the case the simplest scenario would be one where all genus contributions 
gave rise to the same a-dependence. In both case we should be able to test a = 7 + K 
by considering spherical topology. 

In sect. 2 we define the discretized model we are going to use in our attempt to 
verify a = 7 + k. The results of numerical simulations are discussed in sect. 3. 

2 The model 

Let us recall the construction of the sign factor $>(S) for the regular lattice surface 
|3|, 0]. We shall cover the lattice surface by a system of closed non-self-intersecting 
curves, which cross each link ones. This type of covering by curves can be obtained 
by drawing two parallel lines connecting the midpoints of neighbouring links on 
each plaquette of the surface. There are two ways of drawing these lines on each 
plaquette (see fig. 1 in ref. |§) and we have a class of 2 M coverings (M is the total 
number of plaquettes of the surface). Each covering consists of a set of closed non- 
self-intersecting curves Ci,...,C m . Let n(C) be the number of times the curve C 
crosses the lines of self-intersection for the surface S immersed in the cubic lattice. 
Assume we can find a function $(C) such that 

$(C) = (-l) n(c) . (2.1) 

It follows that 

= (-l)« s > = TJ(-1)«(^) = J] $(C fc ) (2.2) 

k k 

where the product is over all curves of the chosen curve system. We can define a 
function $(C) with the above mentioned properties by (this definition makes 
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sense even for a continuum surface): 

$(C) = itr P exp iT 1 ^ . (2.3) 

The matrix R is defined as follows: Consider three orthonormal vectors ef , a — 
1, 2, 3, attached to each point of the curve C : ej is the tangent vector, e\ is normal 
to ej in the tangent plan of the surface and ef is normal to the surface. Then 
R rotates this orthonormal system into a given fixed orthonormal system. If the 
rotation matrix R is in the fundamental representation and denote the Pauli 
matrices this can be written as follows: 

Rr a R~ l = a a , r a = efn (2.4) 

It can be shown that $(C) = 1 if C does not cross a line of self-intersection and -1 if 
it crosses it once. In addition $(C) is invariant under smooth deformations of C since 
it takes discrete values. Only when the deformation of C reaches singular points of 
the immersed surface, where R is not defined, will it be able to change its value. The 
singular points are the possible endpoints of lines of self-intersection of the immersed 
surface and these points are so-called Whitney singularities, which are related to 
7Ti(SO(3)), i.e. elements like ( [Of ). In the case of piecewise linear surfaces (either 
the cubic lattice considered above or the triangulated random surfaces immersed 
in R 3 , to be considered below) we will assume that the curve C is a straight line 
on the flat part of the surfaces. Then the integral (|2.3| ) will only get contributions 
when the curve crosses from one flat piece to a neighbouring one. Let us denote the 
successive flat pieces of surface encountered by the curve by S n , the corresponding 
orthonormal frames by e a (n) and the the rotation matrices related to e a (n) by R(n). 
With this notation ( |2.3|) can then be written as 

«(c)4*n (2.5) 

'l + ef(n)ef(n + l) 



As was shown in ref.0 the sign factor <&(S) does not depend on the specific choice 
of one of 2 M coverings. One can choose the covering such that we have a maximum 
number of closed non-connected curves. These curves will then surround half of the 
lattice vertices separating them into two classes, say "+" and "-" vertices (see fTj). 
Then the sign factor <&(S) can be said to be located at one type of the vertices, say 
the "+" vertices, and by construction we have 

= (-l) n+ (2.6) 
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where n + is the number of Whitney singularities, located at the "+" vertices f] 

Rather than applying the formula (|2.5| ) or ( [2.6|) to the cubic lattice surfaces we 
will assume that there is indeed a continuum surface theory of the 3d Ising model 
at the critical point, and that we are free to use any sensible approximation to 
the continuum surface theory, as long as it admits us to take the correct scaling 
limit. From this point of view it seems reasonable to use triangulated piecewise 
linear surfaces immersed in three dimensions as our statistical ensemble of surfaces. 
Their critical properties coincide with those of the bosonic string in target space 
dimensions d < 1, where the theory can also be solved analytically. We have however 
to choose an appropriate curve system on the piecewise linear surface and for that 
purpose the generic triangulation of a surface is not useful. Rather we consider the 
class of triangulations obtained by gluing squares together as indicated in fig. 1. 
Each square consists of two triangles, one of the diagonals of the square being the 
common link of the two triangles. The four vertices are divided into two groups: 
"-"vertices connected by the diagonal and "+" vertices separated by the diagonal. 
The rule for gluing together links of different squares to form a (triangulated) surface 
is that "+" vertices should be glued to "+" vertices and "-"vertices to "-"vertices 
(see fig. 1). For this class of triangulated surfaces one can now choose a convenient 
curve system which separates the "+" and the "-" vertices. For each square two 
parts of different curves run parallel to the diagonal connecting the "-"vertices. 
Stated differently we can say that the closed curves on the surface surround all 
"+" vertices (see fig. 1). From this point of view the curve system is identical 
to one described above on the cubic lattice, but of course the connectivity of the 
triangulated surfaces are such that they cannot in general be mapped into a cubic 
lattice. The class of triangulations generated this way can be viewed in another 
way: The triangles surrounding "+" vertices form polygons, where the vertices on 
the boundaries are "-"vertices, and precisely one curve from the curve system runs 
around the "+" vertex inside this polygon. The surface is now constructed by gluing 
together these polygons to form a connected closed surface. 

As this class of triangulations is different from the one usually used in discretized 
2d quantum gravity one could worry if they belong to the same universality class, 
which is clearly what we want. One can easily prove this. It is possible to write down 
a matrix model which precisely generates the triangulations we are considering: 

Z(\) = J d<f)d<P*dA e" tr ^^ tr A2+Atr +** A . (2.7) 

4 Note that on piecewise linear surfaces the Whitney singularities have to be located on the 
vertices (the O-simplexes) of the surface. 
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In ( |2.7|) denotes a complex N x N matrix, A an Hermitian N x N matrix and 
to tr cf^ipA we can associate an oriented triangle with one "+" vertex and two "- 
"vertices. The link connecting the "-"vertices corresponds to the A matrix. If 
we first integrate over the Hermitian matrices we reproduce precisely the complex 
matrix model with quartic interactions which is known to have the same ^string 
as the corresponding Hermitian matrix model [K], [|. The integration over the 
Hermitian matrices corresponds creating the squares of fig. 1 by identifying the "-" 
vertices of two triangle. Alternatively the integration over the complex matrices 
will leave us with an Hermitian matrix model with potential a log(l — A). This 
corresponds to gluing together polygons as mentioned aboveQ. 

An immersion in R 3 of a given (abstract) triangulation T of a surface is a 
map from the vertices i e T into G -R 3 , such that links (ij) are mapped into 
straight lines connecting Xi and Xj and triangles (ijk) are identified with the trian- 
gles spanned by Xi,Xj,Xk- This results in a piecewise flat surface S({xi,T}). The 
curve system defined above is by the same procedure mapped into piecewise linear 
curves on S and we can use (|2.2| ) and (|275|) to construct the sign factor Q(S({xi, T}). 
The definition of the model of random surfaces is as follows: 

p -fiN T r 

n») = E ^ttV / 11** 53 (E^) *(%,T}) (2-8) 

TeT u v J ) J ieT i 

where Nj< denotes the number of triangles, C(T) is a symmetry factor for the trian- 
gulation T and T denotes the class of triangulations defined above. In principle we 
can include a summation over different topologies in the sum over triangulations, 
i.e. in the class T, but as discussed above we will here try to test the conjecture 
that surfaces of spherical topology will be important in the scaling limit. This we do 
by comparing the critical exponent a of the 3d Ising model with 7 and k extracted 
from F(p) defined ( |2.8| ), but with T restricted to spherical topology. 

The free bosonic surface theory corresponding to ( |2.8| ) is given by JIT, [T2| , |I3| , |Ti 

FM = E /II dxi 5 3 (E^) e-Xmto->r. (2.9) 

TeT °l J ) J ieT i 

It is in this ensemble we will try to calculate the expectation value ($(S , ))o- We 
expect an exponential fall off like in ( |1 . 1 1| ) , only with Nt instead of the area A. 



5 Note that we have by these arguments shown that to any order in 1/N the m — 2 Hermitian 
matrix model has the same critical behaviour as the m = 2 complex matrix model, a well know 
result which was however not entirely trivial to prove. 
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3 Numerical method and results 



Unfortunately we have at this stage to rely on numerical methods. One such method 
is Monte Carlo simulations. It has been applied extensively to simulations of the 
bosonic string defined by ( |2.9| ). Two types of simulations have been used, a "canon- 
ical" one where the number of triangles Nt are kept fixed, and a "grand canonical" 
one where Nt is allowed to vary|]l4|, [15|, The standard canonical move, the 
so-called "flip" of a link, is not available for us since we consider a restricted class of 
triangulations T. In fact, since we are also interested in a determination of 7 string, 
we have chosen the grand canonical updating scheme, which applies to our class of 
triangulations with a small modification of the standard moves. Each "+" vertex 
is surrounded by "-"vertices and each "-"vertex is surrounded by an alternating 
sequence of "+" and "-"vertices. We have shown the corresponding moves for in- 
serting and deleting "+" and "-" vertices, respectively, in fig. 2. We have checked 
numerically that this algorithm indeed leads to the correct value ^string = — 1/2 in 
the case where the dimension of spacetime is d = 0, but it should be mentioned that 
the extraction of the correct value was more difficult than if we used the standard 
ensemble of triangulations, i.e. all triangulations with spherical topology and length 
of loops of links larger than or equal three. 

Let us now turn to the 3d case. There has been a number of attempts to deter- 
mine 7 string numerically [16, 15, 17, 18, 19]. The conclusion from these simulations 
is that —0.2 < ^string < 0.2. for d = 3. We will not improve these simulations, since 
we will be working with quite small surfaces and the only difference compared to 
these earlier investigations will be that our class of triangulations is different, but 
(at least in zero dimensions) belongs to the same universality class. Our results for 
Istring are compatible with the earlier measurements, indicating that also in three 
dimensions does our class of triangulations give the same results as the unrestricted 
class of triangulations. 

The quantity which has our main interest is (<&n t (S))q, where Nt is the number 
of triangles which constitutes the surface. § Nt (S) is always ±1 but will fluctuate 
wildly, both when we change triangulations and when we update the coordinates Xj, 
and in accordance which our expectations (3>n{S))o falls off exponentially with the 
number of triangles. In practise this means that we have not been able to use surfaces 
with N > 34. These are quite small surfaces, but hopefully large enough^ to reveal 

6 At this point it is worth to recall the history of the determination of ^string by numerical 
methods. Although the first simulations used very small surfaces the qualitative features extracted 
from these simulations are essentially unchanged today, years and thousands of CPU hours later. 
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0.00062 ±0.00033 
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0.00024 ±0.00019 
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0.00040 ±0.00019 
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0.00012 ±0.00020 



Table 1: measurement of {<&n t {S))q 

the qualitative behaviour of the critical exponent k in ($jv(S'))o ~ N K exp(— fiiN). 

In table 1 and fig. 3 we have shown the result of approximately 10 9 sweeps over 
surfaces with iV between 6 and 40. 

The best value one can extract from the data seems to be k — ± 1. The 
determination of k is rather poor and can only be improved by going to larger 
surfaces, which is impossible because of the exponential fall off of the sign factor 
with the size of the surface. We nevertheless find the results encouraging in the 
sense that they are compatible which the relation 7 s t r i ng ±K = a. It would obviously 
be very desirable to be able to find observables which have a better behaviour than 
($Ar)o for large N and which still relate critical exponents of the Ising model to those 
of random surfaces. The situation is the same as for the measurements of Wilson 
loops in lattice gauge theories. The expectation value of the Wilson loops will always 
be exponentially suppressed because we make the importance sampling in the pure 
gauge theory which knows nothing about the static quarks. In the same way we 
here make the importance sampling of pure bosonic string configurations without 
knowledge about the fermionic string which describes the Ising model. If better 
observables could be found it would be possible to go to larger surfaces and thereby 
get more reliable critical exponents and in addition check for the contributions of 
different topologies and orientable versus non-orientable surfaces. 
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Our results indicate that surfaces of spherical topology satisfy the scaling relation 
Istring + k = a near the critical point in a random surface representation of the 3d 
Ising model. We are fully aware that our numerical results are such that "indicate" 
is the appropriate word to use. As discussed above even an improved numerical 
verification of this relation would not prove that only surfaces of spherical topology 
are important in the scaling limit, but it puts emphasis on such a possibility. In 
principle one could dream of settling the problem of dominance of spherical topology 
by numerical simulations since it is possible to perform simulations with surfaces of 
topology different from that of the sphere and also to include non-orientable surfaces. 
However, we need better observables, which do not fall off so fast with the size of 
the surfaces, if we want to be able to make quantitative statements. 

But apart from attempts to improve the numerical simulations we hope that our 
point of view might simulate theoretical attempts to understand better the random 
surface representations of the 3d Ising model. 

Acknowledgement We thank Peter Orland for reading the manuscript and for 
valuable comments. 
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Figure Captions 



Fig.l The matrix model representation of our restricted class of triangulations. Ver- 
tices with open circles represent the "-" vertices defined in the text, while 
vertices with black circles represent the "+" vertices. The dotted diagonals 
symbolize the integration over the Hermitean matrix A^, while the dashed 
line shows the curve C v surrounding the "+" vertex v. 

Fig. 2 The two classes of moves for inserting and deleting "-" vertices (open circles) 
and "+" vertices (black circles). 

Fig. 3 A graphical representation of the data of table 1. We expect for large N: 

log[($ JV )/($ JV _i)] ~ -m + nlog[N/(N - 1)]. 
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